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' Abstract 

■ The Relativistic Hartree Bogoliubov model in coordinate space, with fi- 
nite range pairing interaction, is applied to the description of A-hypernuclei 

I with a large neutron excess. The addition of the A hyperon to Ne isotopes 

■ with neutron halo can shift the neutron drip by stabilizing an otherwise 
unbound core nucleus. The additional binding of the halo neutrons to the 
core originates from the increase in magnitude of the spin-orbit term. Al- 
though the A produces only a fractional change in the central mean-field 

^ . potential, through a purely relativistic effect it increases the spin-orbit 

pH I term which binds the outermost neutrons. 



The production mechanisms, spectroscopy, and decay modes of hypernuclear 
states have been the subject of many theoretical studies. Extensive reviews of 
the experimental and theoretical status of strange- particle nuclear physics can 
be found in Refs. |, H. The most studied hypernuclear system consists of a 
single A particle coupled to the nuclear core. And although strangeness in princi- 
ple can be used as a measure of quark deconfinement in nuclear matter, a single A 
behaves essentially as a distinguishable particle in the nucleus. Theoretical mod- 
els used in studies of hypernuclei extend from nonrelativistic approaches based 
on OBE models for A-N interaction, to the relativistic mean field approximation 
and quark-meson coupling models. However, our knowledge of the A-nucleus in- 
teraction, and of hypernuclear systems in general, is restricted to the valley of 
/?-st ability. 

In recent years the study of the structure of exotic nuclei, produced by ra- 
dioactive nuclear beams, has become one of the most active fields in nuclear 
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physics. The structure of nuclei far from the stabihty hne presents many inter- 
esting phenomena. In particular, in the present work we consider the extremely 
weak binding of the outermost nucleons, large spatial dimensions and the coupling 
between bound states and the particle continuum. By adding either more protons 
or neutrons, the particle drip lines are reached. Nuclei beyond the drip lines are 
unbound with respect to nuclear emission. Exotic nuclei on the neutron-rich side 
are especially important in nuclear astrophysics. They are expected to play an 
important role in nucleosynthesis by neutron capture (r-processes). Knowledge 
of their structure and properties would help the determination of astrophysical 
conditions for the formation of neutron-rich stable isotopes. On the neutron-rich 
side, the drip line has only been reached for very light nuclei. 

It has been recently suggested 0] that a study of A-hypernuclei with a large 
neutron excess could also display interesting phenomena. On one hand such hy- 
pernuclei, corresponding to core nuclei which are unbound or weakly bound, are 
of considerable theoretical interest. Some possibilities for unusual light hypernu- 
clei were already analyzed in the early work of Dalitz and Levi Setti On the 
other hand, one could speculate on the possible role of neutron-rich A-hypernuclei 
in the process of nucleosynthesis. The A particle provides the nuclear core with 
additional binding. Even-core nuclei that are either weakly bound or unbound 
attain normal binding 

In Refs. d^, 1^ we have investigated, in the framework of relativistic mean- field 
theory, light nuclear systems with large neutron excess. For such nuclei the sep- 
aration energy of the last neutrons can become extremely small. The Fermi level 
is found close to the particle continuum, and the lowest particle-hole or particle- 
particle modes couple to the continuum. In Ref. |J7[ the Relativistic Hartree 
Bogoliubov (RHB) model has been applied, in the self-consistent mean-field ap- 
proximation, to the description of the neutron halo in the mass region above the 
s-d shell. As an extension of non-relativistic HFB-theory 0, the RHB theory in 
coordinate space provides a unified description of mean-field and pairing corre- 
lations. Pairing correlations and the coupling to particle continuum states have 
been described by finite range two-body Gogny-type interaction. Finite element 
methods have been used in the coordinate space discretization of the coupled 
system of Dirac-Hartree-Bogoliubov and Klein-Gordon equations. Solutions in 
coordinate space are essential for the correct description of the coupling between 
bound and continuum states. Calculations have been performed for the isotopic 
chains of Ne and C nuclei. Using the NL3 parameter set for the mean-field 
Lagrangian, and DIS parameters for the Gogny interaction, we have found 
evidence for the occurrence of multi-neutron halo in heavier Ne isotopes. We 
have shown that the properties of the lf-2p orbitals near the Fermi level and the 
neutron pairing interaction play a crucial role in the formation of the halo. In the 
present work we essentially repeat the calculations of Ref. [0 for the Ne isotopes, 
but we add a A particle to the system. We are interested in the effects that the 
A particle in its ground state has on the core halo-nucleus. 
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Relativistic mean-field models have been successfully applied in calculations of 



nuclear matter and properties of finite nuclei throughout the periodic table [p!2 
The model describes the nucleus as a system of Dirac nucleons which interact in 
a relativistic covariant manner the isoscalar scalar cr-meson, the isoscalar vector 
w-meson and the isovector vector p-meson. The photon field [A) accounts for the 
electromagnetic interaction. For hypernuclear systems the original model has to 
be extended to the strange particle sector. In particular, the effective Lagrangian 
for A hypernuclei reads 

C = ijj ['^{id - g^u - gppT - eA) -m- g^a) i/j 
-l^,^^'" + \rny - 

+i'A{l{id - gu;AUj) - rriA - gaAcr)iJA. (1) 

The Dirac spinors (■0) and i/ja denote the nucleon and the A particle, respec- 
tively. Coupling constants g^^, g^, gp, and unknown meson masses nio-, m^, nip 
are parameters, adjusted to fit data on nuclear matter and finite nuclei. The 
model includes the nonlinear self-coupling of the cr-field (coupling constants g2, 
93) 

U{a) = ^mla' + ^g^a' + ^gsa^ (2) 

Since the A particle is neutral and isoscalar, it only couples to the a and uj mesons 
(coupling constants g^A and g^jA)- We only consider the A in the Is state, and 
therefore do not include the tensor A - uj interaction.. The generalized single- 
particle hamiltonian of HFB theory contains two average potentials: the self- 
consistent field r which encloses all the long range ph correlations, and a pairing 
field A which sums up the p]5-correlations. In the Hartree approximation for the 
self-consistent mean field, the Relativistic Hartree-Bogoliubov (RHB) equations 
read 

ho-m-X A \fUk\ t:. fUk 



where hn is the single- nucleon Dirac hamiltonian, and m is the nucleon mass. Uk 
and Vk are quasi-particle Dirac spinors, and Ek denote the quasi-particle energies. 
The Dirac equation for the A particle 

[-iaV + (3{mA + g^Acr{r)) + 5'a;A^^°(r)]V'A = CaV'a (4) 

The RHB equations for the nucleons and the Dirac equation for the A are 
solved self-consistently, with potentials determined in the mean-field approxima- 
tion from solutions of Klein-Gordon equations for mesons and Coulomb field: 

[-A + ml] a(r) = -g^ ^ (r)7Vfc(r) - g, a\v) - g, a\v) 

Ek>0 
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-f/.A^l(r)7VA(r) (5) 

[-A + ml] uj%r) = VHr)V,{r) + g^j,i;i{r)Mr) (6) 

Ek>0 

hA + mgp°(r) = Y: y^i^)r3V,{T), (7) 



Ek>0 



-AA^r) = Y VHr)^V,{r). (8) 

Ek>0 

The sums run over all positive energy states. The system of equations is solved 
self-consistently in coordinate space by discretization on the finite element mesh. 
In the coordinate space representation of the pairing field A, the kernel of the 
integral operator is 

A,,(r, r') = ^ E Vabcdir, r')K,,{r, r'). (9) 

c,d 

where V^bcd(r, r') are matrix elements of a general two-body pairing interaction 
and Kcd{r, r'), is the pairing tensor, defined as 

KUr,r'):= ^ f/4(r)y,fc(r'). (10) 

Ek>0 

The integral operator A acts on the wave function Vfc(r): 

(AF,)(r) = E / ciVA,,(r,r')H.(rO. (11) 

In the particle-particle [pp) channel the pairing interaction is approximated by 
the finite range two-body Gogny interaction 



ypp(l,2) = E e"'^''"'"'^^''''' (^* + B^P"" - H^P^ - MiP^P^), (12) 

j=l,2 



with the parameters //j, Wi, Bi, Hi and Mj (z = 1,2). 

As in Ref. the even-even Ne isotopes have been calculated with the 



NL3 [10 1 effective interaction for the mean-field Lagrangian, and the parame- 



ter set DIS [0 has been used for the finite range pairing force. The coupling 



constants for the A particle are from Ref. [13 , where the relativistic mean-field 
theory was used to study characteristics of A, S and H hypernuclei. While the val- 
ues for the g^oY coupling constants were determined from the naive quark model, 
that is g^A = ^dujN', the values of g^y were deduced from the available experimen- 
tal information of hyperon binding in the nuclear medium. For the A hyperon 
(7o-A was fitted to reproduce the binding energy of a A in the Is state of ]^0: 
QaA = 0.62l5(o-Ar. The coupling constant determined from only this experimental 
quantity gives a reasonable description of binding energies in A hypernuclei for a 
wide range of mass number. 
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In Ref. we have shown that the neutron rms radii of the Ne isotopes follow 
the mean-field N^/'^ curve up to N ?5i 22. For larger values of N both neutron and 
matter rms radii display a sharp increase, while the proton radii stay practically 
constant. The sudden increase in neutron rms radii has been interpreted as 
evidence for the formation of a multi-particle halo. The phenomenon was also 
clearly seen in the plot of proton and neutron density distributions. The proton 
density profiles do not change with the number of neutrons, while the neutron 
density distributions display an abrupt change between ^°Ne and ^^Ne. The 
microscopic origin of the neutron halo has been found in a delicate balance of the 
self-consistent mean-field and the pairing field. This is shown in Fig. la, where 
we display the neutron single-particle states If 7/2, 2p3/2 and 2pi/2 in the canonical 
basis, and the Fermi energy as functions of the mass number A. For A<32 (N< 
22) the triplet of states is high in the continuum, and the Fermi level uniformly 
increases toward zero. The triplet approaches zero energy, and a gap is formed 
between these states and all other states in the continuum. The shell structure 
dramatically changes at N> 22. Between A=32 (N = 22) and A=42 (N = 32) the 
Fermi level is practically constant and very close to the continuum. The addition 
of neutrons in this region of the drip does not increase the binding. Only the 
spatial extension of neutron distribution displays an increase. The formation of 
the neutron halo is related to the quasi- degeneracy of the triplet of states If 7/2, 
2p3/2 and 2pi/2- The pairing interaction promotes neutrons from the If 7/2 orbital 
to the 2p levels. Since these levels are so close in energy, the total binding energy 
does not change significantly. Due to their small centrifugal barrier, the 2p3/2 
and 2pi/2 orbitals form the halo. The last bound isotope is ^°Ne. For N > 32 
(A=42) the neutron Fermi level becomes positive, and heavier isotopes are not 
bound any more. 

In the present work we have repeated the calculation of Ne isotopes, but a 
A hyperon has been added to the even-even cores. As one would expect, there 
are no excessive changes in the bulk properties. For example, the neutron rms 
radii are reduced on the average by 2%. Therefore we do not display changes 
in macroscopic quantities, but going to the microscopic level, in Fig. lb we 
illustrate the effect of the A hyperon on the triplet of neutron states that form the 
halo: lf7/2, 2p3/2 and 2pi/2, and on the Fermi level. The energies are displayed 
as function of the core mass number Ag. Due to the extra binding provided 
by the A, the single-neutron energies and the Fermi level are lower. The most 
important effect that we observe, however, is that the Fermi level is negative 
for the isotope ^^^^Ne. Without the A, the nucleus ^^Ne was unbound. The 
presence of the strange baryon stabilizes the otherwise unbound core. This could 
have interesting consequences for the process of nucleosynthesis, especially r- 
processes. It should be noted, however, that heavier nuclei, and in particular 
those with atomic number around Z=40, are important in the rapid neutron- 
capture mechanisms. In a very recent study ||14[, it has been shown that RHB 
theory predicts the last bound Zr isotope to be ^^^Zr (N=98). Therefore we 
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have calculated the ground state of the hypernucleus ^^'^^^Zi. It turns out that 
although the presence of the A lowers the Fermi level, the A-N interaction does 
not have enough strength to stabilize the core Ac=140. But if a single A cannot 
stabilize such a large core, maybe two A particles could. An interesting question 
is whether such objects could be found in the environment in which the processes 
of nucleosynthesis occur. 

It is important to understand the microscopic mechanism through which the 
A binds the additional pair of neutrons in ^^"'"^Ne. In a recent study we have 
used the relativistic Hartree-Bogoliubov model to analyze the isospin dependence 
of the spin-orbit interaction in light neutron-rich nuclei. It has been shown that 
the magnitude of the spin-orbit potential is considerably reduced in drip line 
nuclei. With the increase of the neutron number, the effective one-body spin- 
orbit potential becomes weaker. This results in a reduction of the energy splittings 
between spin-orbit partners. The reduction of the spin-orbit potential is especially 
pronounced in the surface region. In the relativistic mean-field approximation, 
the spin-orbit potential originates from the addition of two large fields: the field 
of the vector mesons (short range repulsion), and the scalar field of the sigma 
meson (intermediate attraction). In the first order approximation, and assuming 
spherical symmetry, the spin orbit term can be written as 

Vs.o. = --^VUr), (13) 
r or 



where Vis is the spin-orbit potential ||T2|, 115 1 



Vis = ^{V-S). (14) 

^eff 

V and S denote the repulsive vector and the attractive scalar potentials, respec- 
tively. rUeff is the effective mass 

meff = m-^{V - S). (15) 

Using the vector and scalar potentials from the self-consistent ground-state 
solutions, we have computed from (|13D - ([15|) the spin-orbit terms for several 
Ne isotopes and corresponding A-hypernuclei. They are displayed in Fig. 2 as 
function of the radial distance from the center of the nucleus. The magnitude 
of the spin-orbit term Vg.o. in Ne nuclei (dashed lines) decreases as we add more 
neutrons, i.e. more units of isospin. The reduction for nuclei close to the neutron 
drip is ~ 40% in the surface region, as compared to values which correspond to 
beta stable nuclei. For the corresponding A-hypernuclei (solid lines) the spin-orbit 
term displays an increase in magnitude of about 10% (smaller as we approach 
the drip line). This is a rather surprising result, as we have seen that the single- 
particle energies and bulk properties seem to be less affected by the presence of 
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the A particle. The effect is purely relativistic and it appears to be strong enough 
to bind an additional pair of neutrons at the drip line. The mean field potential, 
in which the nucleons move, results from the cancelation of two large meson 
potentials: the attractive scalar potential S and the repulsive vector potential V: 
V+S. The spin-orbit potential, on the other hand, arises from the very strong 
anti-nucleon potential V-S. Therefore, while in the presence of the A the changes 
in V and S cancel out in the mean- field potential, they are amplified in Vig. We 
illustrate this effect on the example of ^°Ne and the corresponding hypernucleus 
^^Ne. For the core '^'^Ne the values of the scalar (S) and vector (V) potential in 
the center of the nucleus are -380 MeV and 308 MeV, respectively. For |^Ne the 
corresponding values are: -412 MeV and 336 MeV. The addition of the A particle 
changes the value of the mean-field potential in the center of the nucleus by 4 
MeV, but it changes the anti-nucleon potential by 60 MeV. This is reffected in the 
corresponding spin-orbit term (Fig. 2), which provides more binding for states 
close to the Fermi surface. The additional binding stabihzes the hypernuclear 
core. 

In conclusion, we report results of the first application of the Relativistic 
Hartree Bogoliubov model in coordinate space, with finite range pairing interac- 
tion, to the description of A-hypernuclei with a large neutron excess. In particu- 
lar, we have studied the effects of the A hyperon in its ground state on Ne nuclei 
with neutron halo. Although the inclusion of the A hyperon does not produce 
excessive changes in bulk properties, we find that it can shift the neutron drip 
by stabilizing an otherwise unbound core nucleus at drip line. The microscopic 
mechanism through which additional neutrons are bound to the core originates 
from the increase in magnitude of the spin-orbit term in presence of the A parti- 
cle. We find that the A in its ground state produces only a fractional change in 
the central mean-field potential. On the other hand, through a purely relativistic 
effect, it notably changes the spin-orbit term in the surface region, providing ad- 
ditional binding for the outermost neutrons. Neutron-rich A-hypernuclei might 
have an important role in the process of nucleosynthesis by neutron capture. 

This work has been supported by the Bundcsministerium fiir Bildung und 
Forschung under project 06 TM 875. G. A. L. acknowledges support from DAAD. 
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Figure Captions 



• Fig.l lf-2p single-particle neutron levels in the canonical basis for the Ne 
(a), and A^e + A (b) isotopes. The dotted line denotes the Fermi level. 

• Fig. 2 Radial dependence of the spin-orbit potential in self- consistent solu- 
tions for the ground-states of A^e (a), and A'^e -|- A (b) isotopes. 
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Fig.1 





200 




-300 ' — ' — ■— ' — ' — ' — ' — ' — ' — ' — ' — ' — ' 

2 4 6 

r (fm) 



1 1 1 1 1 1 1 1 1 1 r 




r (fm) 



